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Abstract. In this paper we prove that given a smoothly conformally compact asymptoti- 
cally hyperbolic metric there is a short-time solution to the Ricci flow that remains smoothly 
conformally compact and asymptotically hyperbolic. We adapt recent results of Schniirer, 
Schulze and Simon to prove a stability result for conformally compact Einstein metrics suf- 
ficiently close to the hyperbolic metric. 



1. Introduction 

In 1989, W. X. Shi initiated the study of the Ricci flow on a noncompact manifold by 
proving that there is a short-time solution to the flow starting at a complete metric of bounded 
curvature, and moreover the flow remains in this class. Recently there has been intense 
activity to understand to what extent the Ricci flow preserves other geometric conditions on 
noncompact manifolds, see [2, 6, 14, 12, 13, 21, 27] for examples. In this paper we prove that 
the Ricci flow preserves the set of smoothly conformally compact asymptotically hyperbolic 
metrics in general dimension for a short time. We begin by introducing these metrics. 

Let M n+1 be the interior of a compact manifold with boundary M. Suppose that x is a 
boundary defining function for dM. This is to say that x is a smooth non-negative function 
on M that vanishes to first order precisely at dM. We say that a metric h is smoothly 
conformally compact if h := x 2 h extends to be smooth metric on M. The Poincare ball model 
of hyperbolic space provides an easy example. 

When \dx\j^ = 1 on dM, we may use x to identify a collar neighbourhood of dM in M with 
[0, e) x dM. We then write h as 

dx 2 + h(x) 
x £ 

for a smooth family of metrics h on dM. 

If h is smoothly conformally compact, with \dx\~ = 1 on dM, and (R cc )ijki = huhjk — hikhji 
denotes the curvature 4-tensor of constant sectional curvature +1, then the curvature 4-tensor 
R of h satisfies 



\R + R cc \ h = 0(x), and 



\V { i ) R\ h = 0(x), for all j. 



For this reason conformally compact metrics with \dx\^ = 1 on dM are asymptotically hyper- 
bolic. It is well known that h is complete and of bounded geometry. 
Recall the Ricci flow is the system of equations 

d T g = -2Rcg(r), , , 

9(0) =h. [Ll) 
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As previously mentioned, it follows from [30] that there is a solution to the Ricci flow, g(r), 
with initial metric h for a short time. 

It will be more convenient to study a normalized Ricci flow. Suppose that g N satisfies: 

/ dtg* = -2ngf j -2Rcgf j , 

U» =h. (L2) 

Setting g(x, r) = (l+2nr)g N (x, Jj- log(l + 2nr)) yields a solution to the original Ricci flow. As 
solutions to the Ricci flow and normalized Ricci flow differ by this time rescaling, we see that 
spatial regularity is preserved. Thus a conformally compact and asymptotically hyperbolic 
solution to the normalized Ricci flow yields a conformally compact solution to the Ricci flow, 
with sectional curvatures that depend on time. Moreover it is straightforward to check that 
the conformal infinity is preserved along the flow. 
The first main result of this paper is the following 

Theorem A. If h is smoothly conformally compact and asymptotically hyperbolic then there 
exists a unique smoothly conformally compact and asymptotically hyperbolic solution g(t) to 
(1.2) (and hence a conformally compact solution to (1.1) ^ for a short time. 

The proof of this Theorem proceeds as follows. First we apply the DeTurck trick to obtain 
a system that may be solved by parabolic PDE techniques. Then conditioning the equation 
appropriately we are able to apply a contraction mapping argument to reprove the existence 
(see Theorem 4.7) of a short-time solution to the flow in 0-H61der spaces, which are Holder 
spaces associated to conformally compact metrics that respect the interior geometry. These 
spaces have no tangential regularity at the boundary, so the final part of the argument (see 
Theorem 5.3) proves that the solution is smoothly conformally compact by applying regularity 
techniques modeled on [23]. We have stated our results for smoothly conformally compact 
metrics, but the arguments here extend to certain initial metrics that are polyhomogeneous. 

Given short-time existence for the Ricci flow it is natural to study the stability of the flow 
about fixed points. In general dimension, stability of hyperbolic space under the Ricci flow 
has been studied independently by Li and Yin [20], Schniirer, Schulze and Simon [29], and 
Bamler [6]. The second main result of this paper extends the stability result of Schniirer, 
Schulze and Simon to certain small Einstein perturbations of the hyperbolic metric. By an 
77-admissible Einstein metric we mean an at least C 3,a conformally compact Einstein metric 
h on B n+1 that satisfies both a global curvature bound 

sup \R + R cc \ h < 7], 

and that the Yamabe invariant of its conformal infinity is positive. Note the existence of such 
metrics follows from work of Graham and Lee [10]. In particular any smooth Riemannian 
metric h on S n sufficiently close to the round metric in an appropriate C k ' a norm is the 
conformal infinity of a conformally compact Einstein metric h on the unit ball. 

Before stating our second main result, we introduce the function spaces defined in [29]. For 
an interval I C [0, oo), let M k (M n+l , I) (resp. M k oc (M n+1 , 1)) denote the space of sections g(t), 
t E / of metrics on B n+1 which are C k (resp. C k oc ) on B n+1 xl, with covariant derivatives taken 
with respect to h. The space A4g°(l n+1 , 1) will denote metrics in M°{M n+1 , 1)nM^ 0C {M n+1 , 1) 
which are smooth for positive times and uniformly bounded in C k when restricted to time 
intervals of the form [5, oo), 5 > 0. 

Adapting the work of Schniirer, Schulze and Simon we prove 
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Theorem B. Let n > 3. There exists an rj(n) > so that for any n- admissible Einstein 
metric with < r) < rj{n) the following holds. For all K > there exists e\ = e\{n,K) > 
where if go is a ftA° metric close to h in the sense that 

/ \g - h\ 2 h dvol h < K, 

and 

sup \g - h\ h < ei. 

1"+! 

Then there exists a long-time solution g(t) G M X) (M n+1 , [0, oo)) to the normalized Ricci- 
DeTurck flow (with initial metric go) such that 

sup \g(t) - h\ h < C(n,Er)e _I (^*. 

Moreover, g(t) — > h exponentially in C k as t — > oo, for all k E N. 

In the above theorem, unlike in [29], we have transcribed the dimension to n+1 to match the 
convention of the rest of the paper. Note also that in view of the main regularity result of this 
paper, if go is smoothly conformally compact, then g{t) remains smoothly conformally compact 
for finite time. The limiting conformally compact Einstein metric need not be smoothly 
conformally compact. 

We take this opportunity to mention two related papers. First, recent work by Hu, Qing 
and Shi [12] proves the Ricci flow preserves a certain class of asymptotically hyperbolic metrics 
for a short-time. These metrics are defined by curvature decay conditions and, as shown in [4] 
and [12], are conformally compact of only a limited regularity. Hu, Qing and Shi subsequently 
prove an interesting rigidity result. On the other hand, in view of the applications of smoothly 
conformally compact metrics to geometry and physics (see for example [7] and references 
therein), it is natural to study the Ricci flow in the smooth conformally compact setting. 
Second, the author and Helliwell have recently proved short-time existence results for higher- 
order geometric flows on compact manifolds [5] . We observed that many short-time existence 
results depend only on the special algebraic structure of the flow. Both [5] and the present 
paper were developed in parallel, and were inspired by recent work of Koch and Lamm [15]. 
The short-time existence of the Ricci flow we give here, while in the setting of conformally 
compact metrics, may be regarded as a concrete application of the ideas in [5]. 

After the present paper was accepted for publication Qing, Shi and Wu [28] posted a 
preprint that also studies the Ricci flow of conformally compact Einstein metrics. Their 
technique differs from the one presented here in that they use maximum principles to obtain 
regularity along the flow. They are also able to prove interesting perturbation results for 
conformally compact Einstein metrics. 

This paper is structured as follows. In Section 2, we outline the DeTurck trick and reduction 
of the flow to a parabolic system. In Section 3, we define function spaces and outline the main 
results from linear parabolic theory on conformally compact manifolds. This theory is based 
on the edge and heat calculus for 0-operators that appears in [22] and [1]. In order to not 
distract from the main Ricci flow argument, we have kept this section short and instead 
sketched several of the proofs of the analytic results in the Appendix. In Section 4, we 
condition the flow equations and provide the contraction mapping argument. We discuss the 
regularity argument in Section 5, and the stability argument in Section 6. Finally, in the 
Appendix we provide sketches for the various analytic facts quoted in Section 3. 
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2. Preliminaries 

As is well known, the Ricci flow is not a parabolic system due to the diffeomorphism 
invariance of the Ricci tensor. We will break this invariance using the standard DeTurck trick. 
Choosing the initial smoothly conformally compact and asymptotically hyperbolic metric h 
as the background metric, and writing all Christoffel symbols and curvature quantities with 
respect this metric with tildes, we define a time dependent vector field 

w k = <T (r k pq - vj q 

The normalized Ricci-DeTurck flow is given by 

d t gij = -2ngij - 2Rc g^ + ViWj + VjW t 



g(P)ij =hij. (2 - 1} 

Standard computations, for example given in [30], show that this flow may be written 
dtgij - g ab V a V b9ij + 2n 9ij + g ab g ip hPm jaqb + g ah g jp h™R iaqb 



-\g ab g pq \ Vi9paSjgqb + 2V a g jp V q g ib - 2V a g jp V b g iq 

~ ~ ~ ~ \ (2.2) 

- 2Vjg pa V b g iq — 2Vig pa V b gj q 



g(0) =h. 

From this equation we see the Ricci-DeTurck flow is a quasilinear parabolic system for the 
metric. 

Once we prove short-time existence of a smoothly conformally compact solution g to the 
Ricci-DeTurck flow, the time-dependent vector field W k will have coefficients smooth up to 
the boundary of M and vanishing to first order there. If (f>t denotes the flow generated by 
W k , then g = <\>* t g is a solution to the normalized Ricci flow [9]. It is straightforward to see 
that g is smoothly conformally compact and asymptotically hyperbolic. 

Finally, we only prove the existence of a short-time solution to the Ricci flow. The unique- 
ness assertion in Theorem A follows from the work of Chen and Zhu [8] . 



3. Parabolic theory on conformally compact spaces 

In this section we outline linear parabolic theory for uniformly degenerate operators on 
conformally compact manifolds. We just state the results we need here; sketches of proofs are 
deferred to the appendix. The primary references for the material in this section are [22] and 

[!]• 

Let (M, h) be a smoothly conformally compact asymptotically hyperbolic manifold as de- 
fined in the introduction. Suppose that x is a boundary defining function and that {y 1 , • • • , y n } 
are coordinates on the boundary, extended to be constant in x. We will refer to these coor- 
dinates as background coordinates. Note that in order to avoid additional notation we will 
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occasionally use the convention x° = x,x l = y l ,i = 1 • • • n to generically refer to the coordi- 
nates as x l when we do not need to keep track of normal/tangential directions. The metric h 
decomposes as 

_ dx 2 + h ab (x,y)dy a dy b 

X A 

where the components of h are smooth up to the boundary. 
The ^-vector fields are generated by 

{xd x ,xd y i,- ■ ■ ,xd y n] , 

and form the basis of a vector bundle, the 0-tangent bundle °TM. We will also have occasion 
to discuss b-vector fields, which are generated by 

{xd x ,d y i,--- ,d y n} . 

We will denote the space of smooth functions on M by C°°{M) and functions smooth up 
to the boundary by C°°(M). The vector bundle of symmetric 2-tensors on M will be denoted 

S 2 (M). We will use ^ and — as the preferred basis for this bundle. 

An operator L on functions is uniformly degenerate of order m if in local coordinates it is 
given by: 

L= ^2 a jA x ^y^ t )( xd ^y( xd v)' 3 - 

j+\P\<m 

where the coefficients aj t p are at least continuous up to the boundary. In order to use Albin's 
heat calculus, we require that aj^ be smooth up to the boundary and independent of time. 

The principal symbol of a uniformly degenerate operator L is a homogeneous polynomial 
on °T*M given by 

°a(m,n)= £ a h ^rf. 

j+\fS\=m 

We say that L is elliptic if °cr(L) is invertible away from (£, rj) = 0. 

For the Ricci flow analysis, we will have to deal with systems of equations as our operators 
will act on the vector bundle of symmetric two tensors. An operator between tensor bundles 
E and F is uniformly degenerate if in local coordinates it may be written as a system: 

{Lu)i= ^2 ( a j^)i( xd x) 3 (xdy)^u k . 
j+\P\<m 

where the coefficients aj t p are now entries of a dim F x dim E matrix that is at least continuous 
up to the boundary. The principal symbol is defined as before. We will not need to consider 
the most general notions of ellipticity for systems as the Ricci flow system (2.2) is 'diagonal' 
at top order, i.e. (aj,p)\ = (o>j,p) ■ 5 l k . From this we can see that all coupling occurs at lower 
order. We now say that L is elliptic if dim F = dim E and the symbol is invertible away from 

= o. 

For the remainder of this section we suppose that L is a second order uniformly degenerate 
elliptic operator with diagonal principal symbol. 

3.1. Function spaces. We work in the 0-H61der spaces defined for example in [19, 22, 24]. 
We describe the anisotropic version of these Holder function spaces, and refer the reader to 
the references for the purely spatial version. For any manifold M, the notation Mt will denote 
the cylinder M x [0, T). Fix a smoothly conformally compact metric h, which in the Ricci flow 
analysis, will be the initial metric. We assume a covering of M by background coordinates 
has been fixed. 
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Cover M by a Whitney decomposition of countably many uniformly locally finite coordinate 
balls Bi with centre (xi,yi) and radius \xi. We will consider the product of each ball with a 
time interval [0, T). For any < a < 1, consider the norm 

(x + x') a \u(x, y, t) — u(x', y', t)\ 

I M \a,% ■= Moo + SUp < SUp : — — — 

\u(x,y,t) - u(x,y,t')\ \ 

| SLIP , / . 

{x >y ,t)^(x,y,t')e(Bi) T \t ~ t'\ a / 2 J 

The prefactor x + x' comes from using the euclidean metric in background coordinates instead 
of the intrinsic ^-distance, see [24]. Note that we may also use an affine map ipi : Bt —> (-Bi)t 
from a fixed standard cylinder Bt to define these norms. 

Let C e ' 2 (Mr) be the closure of C°°(Mt) with respect to this norm. We define 

C e ' 2 (Mr) to consist of all functions u such that {dt) 1 {xd x ) 3 {xd y f u G C e ' 2 (Mr) for 

all 2i + j + |/3| < A;. Note that unlike C°°{Mt), the spaces Ce +a ' 2 (M T ) and even 
^'"(Mr) have no tangential regularity at the boundary. We also weight these spaces: 

k+a k+a k+a k+a 

u G x v C e ' 2 (Mt) if and only if u = x u v for some v G C e 2 (Mr). 

We will also need Holder spaces of tensors. As previously stated, we use the vector fields 
xd x and xd y b and covector fields dx/x and dy b /x as a basis for bundles of tensors, and with 
this convention V h involves only derivatives by the 0-vector fields. In this way a section of a 
tensor bundle is an element of a Holder space if and only if its components are. Furthermore, 
for j < k 

I,. k-\-a a_ 

(V h ) j : x v C e ' 2 (Mr; E) — ► x v C e ' 2 {M T - E ® °T*M). 

Finally, in what follows since we always deal with the bundle of symmetric 2-tensors, we will 
not explicitly mention it in the notation. 

In our final regularity argument we will need weighted Holder spaces that allow for tangen- 
tial regularity. Following [25] we introduce a scale of spaces C fc+a ' 2 >'(Mr) which consists 

of elements u G C e 2 (Mr) such that d s y u G C e ' 2 (M T ) for < s < /. This is to 
say that up to I of the xd y 0-derivatives may be replaced by the tangential d y 6-derivatives. 

We weight these spaces as before. Note that C k+a ^'°{M T ) = d +a '^~ (M T ), the 0-H61der 

space, and that C k+a > 2 > k (M T ) is a Holder space of k 6-derivatives. 

In [19, 22], elliptic estimates in 0-H61der spaces are proved from scaling and classical interior 
elliptic estimates on the balls Bi, as the pullback of a uniformly degenerate elliptic operator 
under ipi becomes uniformly elliptic. Similarly we may obtain parabolic estimates from scaling 
and classical parabolic estimates. In particular we have the following regularity result, see 
[16, Theorem 8.11.1, Theorem 8.12.1] for the classical parabolic statements. 



Proposition 3.1 (Parabolic regularity). Let L be a second order uniformly degenerate elliptic 
operator. Suppose that D~<a jtl3 G (M T ) for | 7 | < k, and f G Ce^{M T ), U 1 ^ G C*{M) 



2+a +a 

for all I7I < k. If u G C e 2 (Mr) is a solution to (d t — L)u(£,i) = f{C,t) then D^u G 
d +a ^{M T ) for all | 7 | < k. 
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3.2. Parabolic Schauder estimates. We now state the main facts from linear parabolic 
PDE theory that we need. We will be interested in the following problem 

(d t -L)u(C,t) =f((,t) 

u(C,0) =0, {6A > 

The basic result is 

Theorem 3.2. Suppose L is a second order uniformly degenerate elliptic operator with 

a - 

time-independent coefficients. For every f G x^Ce 2 (Mt) there is a solution u to (3.1) in 



J2+a 



■2 + it 



x^Ce ' 2 {Mj^). Moreover, u satisfies the parabolic Schauder estimate 

\\u\\ 2+a ,2+. <K\\f\\ a (3.2) 

The Schauder constant K that appears in the statement depends on T but remains bounded 
as T — > 0. Please see the appendix for a proof. 

3.3. Mapping properties of the heat operator. Given the homogeneous Cauchy problem 

(d t -L)u(C,t) =0 

u(C,o) = <MC), [ ' 

let A denote the heat operator such that takes <p to the solution of this problem, i.e. 
(A<f))(C,t) = u(C,t). We also use the notation that A = e tL . In the appendix we describe 
how A is given by an integration against a specific polyhomogeneous distribution on a certain 
manifold with corners that covers M x M x R + . The mapping properties of such operators 
follow from the asymptotics at each of the boundary hypersurfaces. A key result that we will 
need is that if Vf, is a b- vector field and A is a heat operator, then the commutator [A, V&] has 
the same asymptotics as A, and will enjoy the same mapping properties. See Proposition A. 6 
for a precise formulation. 

Let H denote the following time convolution of the heat operator 

(Hf)((,t)= [ t e^ L f(-, S )ds. 
Jo 

This operator provides a solution to the inhomogeneous Cauchy problem with zero initial 
data. The precise mapping properties we need are given in the following 

Proposition 3.3. (see Corollary A. 4) If e x^C°°(M) and f G x^C°°(M T ) then 

(1) A<j> G x^C°°(M T ). 

(2) Hf G x^C°°(M T ). 

Once again we defer the proof to the appendix. 

4. Short-time existence 

In this section we prove short-time existence of a solution to (2.2) in the 0-H61der spaces. 
This is based on a contraction mapping argument. 

We begin by making several observations that will be needed later. Let E = R + R cc 
be the curvature 'error' tensor for the conformally compact metric h, where R cc denotes the 
+ 1 constant curvature 4-tensor. By our convention for function spaces, if h is smoothly 
conformally compact then 

h = h L^^ec^iM), 

XX 
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where we recall the convention that (x°, x l , • • • , x n ) = (x, y 1 , ■ ■ ■ , y n ). 
We also have E G xC°°(M;T 4 M) and thus E G xC™(M;T 4 M). 

k+a k+a 

We need an expansion for the inverse of the metric. Suppose that v G xC e ' 2 (My) with 
sufficiently small norm, then the symmetric 2-tensor h + v will be invertible and (h + v ) _1 G 

k+a fc + a 

C e ' 2 (Mt). Furthermore, we document a useful expansion 

(h + v f b = h ab - h al h bm v ml + (h + v) bl h am hP q v lp v mq . (4.1) 

4.1. Conditioning the Ricci-DeTurck system. Here we pursue short-time existence of 
the normalized Ricci-DeTurck flow. We will look for a solution of the form 

9ij(x,y,t) = hij(x,y) + v ij (x,y,t) 

k+a k+a 

where V{j G xC e ' 2 {Mt)- The system (2.2) for v may be written in the following way, 
which will facilitate treating the quasilinear system with a contraction mapping argument. 
Here we handle the quasilinearity as a quadratic error. 

= d t Vij - h ab V a V b Vij -((h + v) ab - h ab ) V a V bVl] +Mh + v)ij 

-(h + v) ab (h + v) ip hP0R jaib -(h + v) ab (h + v) jp h™R iaqb (4 2) 

+[(h + v)- 1 *(h + v)" 1 *Vv* Vv]ij, 
v(0) = 0. 

Note that in this expression we have switched curvature sign conventions from [30] . Shi lowers 
an index in the curvature tensor to the third slot whereas I lower to the fourth slot. The 
asterisk denotes linear contractions whose precise formula is unimportant for what follows. 
Let us introduce notation for some of the terms above. Define 

(T\v)ij := (fh + v) ab - h ab ) V a V fe% -, 

(T 2 v)ij := 2n(h l3 + % ) + (-(h + v) ab (h + v) ip h™R jaqb - (h + v) ab {h + v) jp h™R iaqb ^ , 
(T 3 v)ij := (h + v)" 1 *(h + v)" 1 *Vv* Vv. 

We begin by studying the various mapping properties of the terms of this equation. Much 
of the argument depends on the special algebraic structure of these equations. We introduce 
the following notation. We will say various terms are Q(v) if they are linear combinations of 
contractions of bounded tensors with either v or its first two /i-covariant derivatives. We will 
loosely refer to this dependence as being 'quadratic', and we will make precise the estimates 
we need at the end of this section. Note that indices on Q index the term of origin in the 
decomposition above. 

Lemma 4.1. 

T\v = Qi(v), and 
Qi : x v Ce 2 (M r ) — > x 2v C e + ' 2 (M T ). 
Proof. We begin by applying the expansion for the inverse in equation 4.1 
(TitOy := ({h + v) ab - h ab ) V a V bVij 

= (h*b _ h *l h b™ Vml + ( h + V fh am h^V lp V mq - h ab ) V a V b Vij 

= h^ 1 * h^ 1 * v * V 2 v + (h + v)" 1 * /i" 1 * hT 1 * v * v * V 2 v, 
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which shows the expression is quadratic in v. Noting that h~ l G C£°(M) and v G 

fc I q- fc+a 

Ce' 2 (Mr), we see that while we lose two O-derivatives we gain decay in x, i.e. 

Qiv G x 2v Cl 2+ "' 2 (M T ). □ 

The expression for T 2 simplifies considerably. Note that in order to recognize the Lich- 
nerowicz Laplacian below we will not surpress indices in the expression that follows. 

Lemma 4.2. 

(T 2 v)ij = 2Eij + 2nvij + VipRc/ + Vj P Rc i P + 2v ml R m ij 1 + Q 2 (v)ij. 

1. I „■ k+a j i k + a 

Q 2 : x u Ce ' 2 (M T ) — > x 2 ^C e fc+ ' 2 (M T ). 

Proof. By applying the expansion for the inverse to terms in T 2 we find the expression contains 
inhomogeneous terms as well as terms linear in v which we must separate from the main 
expression. In particular, considering one of the constituent terms in T 2 we find 

-{h + v) a \h + v) ip h™R jaqb = -{h + v) a \h + v) ip R ja \ 

h ab - h al h bm v ml + (h + v) bl h am h™v lp v mq ) (hi P + v ip )R ja p b 
= Rcij + Vi p Rcj P — v m iR m ij 1 + [h^ 1 * h^ 1 *v*v* R]ij, 
where R in this calculation denotes the (1, 3) tensor. One may check that the final quadratic 

h I a- hJra fc 1 q- fe ~*~ a 

contraction terms map x u C e ' 2 (Mr) — > x C e ' 2 (M^). 
Consequently, 

-{h + v) ab (h + v) ip h™R jaqb -(h + v) ab (h + v) jp hPiR iaqb 

= 2R Cij + VipRc/ + %,#c/ - 2v m iR m ij 1 + Q 2 (^)y 

Note that by the curvature asymptotics Rc^j = —nhij + where E^ G xC°°(M). There- 
fore, re-assembling T 2 we find 



(T 2 v)ij = 2E i:j + 2nvij + v ip Rc/ + VjpRc/ - 2v ml R m ij 1 + Q 2 (v 



□ 



The third term requires no additional conditioning. 
Lemma 4.3. 

(T 3 v)ij = Q 3 (v), and 

1 I n k-\-a j -I 1 fc— 1-j-q 

Q 3 : x"CT* 2 (M T ) — > x 2l/ C7 e ' 2 (Mr). 

The preceding lemmas allows us to condition the equation for v further. We now move 
the terms linear in v to the other side of the equation. We also have from [9] that the term 
h ab Va^bVij is the rough Laplacian on 2-tensors. In fact, we see the linear elliptic part of the 
equation is the Lichnerowicz Laplacian on 2-tensors, 

L = A^Vij + 2n Vij = h ab V a V bVi j + v ip Rc/ + vjpRc/ - 2v m iR m ij 1 + 2n Vij . 



We may write: 



%(0) =0. 
For the remainder of the argument we drop indices 



d t Vij - (Lv)ij = Qvij + 2Eij, ^ ^ 
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To summarize the argument so far, we have conditioned the flow equations to recognize a 
strongly parabolic equation for the metric. As the quadratic terms Q depend on v and up to 
its first two covariant derivatives in a polynomial fashion, there is a constant C > depending 
on the algebraic structure of Q such that for all u,v G x' l C 2+a, 2(MT), 



Q(v)\\ a % <C\\v\\ 2 2+a a , (4.4) 
x»C7^{M T ) 



< C max < Hull 2 + a a , I \v 1 1 2+ a a >||u — vll 2 +a a . (4.5) 

Note in these estimates that we are relaxing control of one time derivative. This will 
facilitate the contraction mapping argument given in the next section. Note also that this 
part of the argument will not explicitly use the gain of decay by Q. 

In the regularity argument of Section 5, we will need this additional decay. We conclude 
this section with the following lemma. 



Lemma 4.4. All of the quadratic mapping terms satisfy 
Moreover, we have 



lc4-n fe + a r, b — O-l-n fc — 2 + q 

Q : x u Ce 2 (M r ) — »• x 2v C e ' 2 (Mi 



k4-n k+a I hA-n k+a 

if w = w ' + w", where w' G x v C e ' 2 ' (Mr) and w" G x^C e ' 2 (M T ), (u < fj,) 



then 



L o i , k — 2-\-a 7_o h, o i n k — 2-\-a 

Q(w)£x 2u Ce ' 2 ' (M T ) + x' ,+v Cr t ' 2 (M T ) 



k+a k+a I 

• if w = w ' + w " > where w' G x u C°°(M T ) and w" G x^C e ' 2 ' (M T ), (y < fi) then 
Q(w) G x 2v C°°{M T ) + x» +v d~ 2+a ' h= ^' l ~ 2 (M T ). 

Proof. The first mapping property stated follows from the previous lemmas. We need only 
check the final mapping properties. These are straightforward to check as V acts by 0- 
derivatives thus preserves the order of decay. The explicit contraction structure of each of the 
terms that form Q are: 

Qiv = h^ 1 * h^ 1 * v * V 2 v 

Q2V = h^ 1 * h^ 1 * v * v * R 

Q3V = (h + v)^ 1 * (h + v)" 1 * Vv * Vv. 

Inserting w = w' + w" into the expression we find the cross terms have the decay expected of 
w' *w". □ 

4.2. The contraction mapping argument. We now explain the contraction mapping ar- 
gument that leads to short-time existence for equation (4.3). Write the heat operator for 
dt — L as e tL . Apply Duhamel's principle to (4.3) to get an equivalent integral equation 

v(t) = C e^ L {E + Q{v)) ds . (4.6) 







Note the definition of the map ^ in the displayed equation above. 
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2+a — 

For a parameter \i and T to be specified, define a subspace -2^,t of xC e 2 (-Mt) by 
Z„, T = \uexC 2 e +a ^(M T ) :u(x,0) =0,\\u\\ ^+-.§,^ < M-} ■ 

[ xC e J (M T ) J 

This is a closed subset of a Banach space. 

Suppose that u £ Z^ t T, it follows that v = is a solution to 

(d t - L)v = Q(u) + E, 
v(0) = 0. 

As Q(u) + E £ xC a '%(M T ), the Schauder estimate implies v £ xCl +a,1+1 (M T ) C 

2+a - 

xC e ' 2 (Mr), and so 

* : Z^ T ^ xCe +a 'HM T ). 

We now prove that \& in fact maps -2^,t to itself and is a contraction for /j, and T sufficiently 
small. 

Lemma 4.5. ^ : Zn t T — > Z^t for fi and T sufficiently small. 
Proof. To begin, let u £ Zu,T and set 



V! := / e^-'^Q^ds 
J o 



u 2 := 

Consider ux. This is a solution to 



e {t - s)L Eds. 



{d t -L) Vl =Q(u), 
vi(0) = 0. 

The Schauder estimate, followed by the estimates for Q given by equations (4.4) and (4.5) 
gives 

|W U 2+ -v T ) - lw U 2+ - i+ v T ) 

11 xC e (Mx) 

< KC\\U\\ 2 2 + a a 

< KCfl\\ U \\ n 2+a,% . 

Taking /i sufficiently small allows us to force KCfi < |. So ||t>i|| ^ 2 + a ,§^ ^ < ^- Note that 

this same \i works if we shrink T. 

Regarding i>2, note that this is a solution to 

' (d t -L)v 2 =E, 
v 2 (0) = 0. 

We recall that E and the coefficients of L are smooth, time-independent and have bounded 
0-derivatives of all orders, and so by parabolic regularity any finite number of derivatives of 
V2 are bounded. Fixing any £ we may write 



V2((,t) 



f E(() + Lv 2 ((,s)ds, te [0,T). 
Jo 
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2+a 1+ — 

We may now estimate the xC e ' 2 (Mt) norm of V2- The L°° norm of spatial derivatives 

may be controlled through the Schauder estimates by the norm of E, and can be made as small 

as we like by choosing T sufficiently small. Further, as the time derivative of V2 is bounded 

a - 

and V2{x, 0) = 0, the C e ' 2 (M-jp) norm of V2 can be made arbitrarily small by choosing T 
sufficiently small. We conclude for T small enough 

ii ii ^ ^ 

W V2 \\xC 2+a '%(M T ) ~ 2' 

Thus : Z^ T — ► Z^ T for t £ [0, T\. □ 

Lemma 4.6. For the fi and T specified in the previous lemma, ^ : Z^t — > Z^ t T is a 
contraction. 



Proof. Schauder's estimate applied to — implies 

\\^fU — ^>V\\ 2 + a St 

< U^U - tyv\\ 2+a,l+& 

< K\\Qu-Qv\\ 

<KCn\\u-v\\ xcl+a ,« (MT) . 

Where K and C are the same constants from the previous proof. Consequently KC[i < ^, 
and V& is a contraction. □ 

We are now ready to prove the existence of a solution to the Ricci-DeTurck flow with full 
0-regularity. 

Theorem 4.7. If h is a smoothly conformally compact metric, then there exists T > and a 
solution g £ CT ,C °{M T ) to (2.2). 

Proof. The existence of a solution to (4.3) in Z„t follows from the Banach fixed point theorem. 
The Schauder estimate applied to the fixed point equation shows that the solution lies in 

2+a 2+a 

C e ' 2 {Mt). This short-time solution yields a solution in the same space to the Ricci- 
DeTurck flow by taking g = h + v. We now improve the regularity by using a bootstrap 
procedure, applied to the system (2.2). We may write this abstractly as 

2 

dt9+ a ^9) D ^9, 

|/}|=o 

where the coefficients ap at worst satisfy D^ap £ C a, 2 (Mt), for I7I = 1. By parabolic 

regularity (c.f. Proposition 3.1) we conclude D~> g £ C e 2 (Mr) for all I7I = 1, which 
allows us to improve the spatial regularity. By bootstrapping, and then using the equation to 
improve regularity in time, we find g £ C^ ' 00 (Mt). □ 

5. Beyond 0-regularity 

In the previous section we proved short-time existence of the Ricci-DeTurck flow starting 
at a smoothly conformally compact metric. The solution was constructed in 0-H61der spaces 
and is smooth in time and 0-derivatives. We can expect more as the inhomogeneous terms in 
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equation (4.6) are smooth up to the boundary with respect to background derivatives. In this 
section we prove the solution remains smoothly conformally compact on the entire interval of 
existence. The arguments of this section are modeled on the arguments in [23]. 
We begin by writing (4.6) more compactly as 

v = HE + HQv, (5.1) 

where H is the time convolution of the heat operator appearing in (4.6). 

We begin by documenting the mapping properties of H on the C k+a ' z ''(My) spaces. 

Lemma 5.1. 

H : x^C k+a '^' l {M T ) — > x^C k+a ^\M T ),l < k. 

Proof. When I = 0, the mapping properties follow from the Schauder estimate, Theorem 3.2. 
The new content here is that H preserves tangential regularity. 

Suppose that / G x^C k+a '^2~' 1 (Mj')- Taking an arbitrary 6-derivative of Hf, we write 

dy(Hf)=H(d y f) + [H,dy]f. 

Since d y f G x"C k - 1+a > h= * t *(MT), and H : x^C k - 1+a '^ (M T ) — >• x'*c£ +1+a, * ± = ±2 (Mr), 

k+a k+a 

the first term lies again in x^C e ' 2 (Mt)- For the second term, by Proposition A. 6, 

fe4- a k + a 

[H, d y ] has the same mapping properties as H, and again maps x^C e ' 2 (Mt) to a 
subset of xt*Ce +a,k ^(M T ). This implies that d y (H f) G x^Ce +a '^ (M T ), and so if/ G 

h4-n k+a 1 

C e + ' 2 ' (Af T ). 

The remainder of the proof proceeds by iteration. □ 

The next proposition shows that a solution to the normalized Ricci-DeTurck flow is smooth 
in ^-derivatives. 



k+a 



k+a 



Proposition 5.2. Let g be a solution to (2.2) in C e 2 (M?), with h smoothly conformally 
compact. Then g(t) lies in C k + a '^> k for all t G [0, T). 

Proof. Consider the term HE in equation (5.1). As h is smoothly conformally compact, E G 
xC°°(My). Now H preserves polyhomogeneity via Proposition 3.3, and so HE G xC°°(Mt). 
Thus we need only focus on the second term. 

In order to handle the term HQv we take advantage of the improved decay of Qv. If 

v G xCe 2 (M T ), then Qv G x 2 C e 2 (Mr), which H maps to x 2 C e 2 (Mr). 

Consequently, taking any 6-derivative of HQv yields 

dyHQv = x" 1 [xdy(HQv)) G xC e ' 2 (Mr), 

which shows a gain of one tangential derivative so that 

HQv G xC k+a ^'\M T ). 

This argument iterates fc-times and we indicate the next step. Apply equation (5.1) and 
Lemma 4.4 to obtain 

v = HE + HQ(HE + HQv) G xC°°{M T ) + x 3 C k+a ' ! ^' 1 (M T ), 
and so v G xC k+a '^> 2 {M T ). □ 
We now state the main regularity result of this paper. 
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Theorem 5.3. Let g be a solution to (2.2) in C^°' qo (Mt), with h smoothly conformally 
compact. Then g{t) is smoothly conformally compact for all t G [0,T). 

Proof. By Proposition 5.2 we have that g(t) is fully tangentially regular, i.e. v G 
xC° ' 00, °°(Mt)- It remains to show that v is smooth up to the boundary, i.e. v G xC°°(Mt). 

We will now use the structure of the heat kernel as a polyhomogeneous distribution to 
prove polyhomogeneity of v. Given that v G xC 00 ' 0C ' 00 (Mt) satisfies 

v = HE + HQv, 

We see that Qv G x 2 C°°'°°'° (Mt), and so we may decompose v as 

v = v > + v " e xC°°(M T ) + x 2 C co ' 00 ' 00 (M T ). 
We now insert this back into (5.1). Using Lemma 4.4, we find 

Q( v ' + v ") = X 2 C°°(M T ) + ^ c oo,oo,oo (Mt ^ 

Equation (5.1) now lets us conclude 

v G xC7°°(M T ) + x 3 C 00 ' 00 ' 00 (M T ). 

Iterating we conclude v G xC°°{Mt). 

Finally, we have proved v G xC°°(Mt), i.e. that 

dx l dxi 

V = XVij , 

xx 

where Uij is smooth up to the boundary, and we remind the reader of the convention 
(x°,x , ••• ,x n ) = (x, y , ••• ,y n ). So now x 2 v = xv i jdx' L dx^ and consequently g = h + v 
is smoothly conformally compact. □ 

We note that g = x 2 g = {hij + xvij)dx l dx^ , so that \dx\j = 1 on dM, so that g is 
asymptotically hyperbolic. Combining Theorem 4.7 and Theorem 5.3 completes the proof of 
Theorem A. 

We conclude this section by remarking that it is possible to relax the condition that the 
initial metric h be smoothly conformally compact. Our entire argument applies to certain 
polyhomogeneous initial metrics as well. The heat kernel analysis from the Appendix and the 
heat kernel mapping properties extend to this case in a straightforward manner. 

6. Stability about admissible Einstein metrics 

In this section we adapt the arguments from [29] to prove stability of the normalized Ricci- 
DeTurck flow near 77-admissible Einstein metrics. 

The main idea is to replace the background hyperbolic metric on the ball used in [29] with 
an 77-admissible Einstein metric, h. The existence arguments in [29, Theorem 2.4] works with 
a such a complete Einstein metric in place of the hyperbolic metric. In order to obtain the 
long-time existence and convergence result, we must adapt the L 2 estimate for the difference 
of the solution of the flow with the background metric due to the 'curvature error' terms, 
whose magnitude is measured by 77, that arise. There is sufficient room in the main estimate 
to allow for metrics with n sufficiently small. To this end, we replace [29, Lemma 2.2] with 
the following 
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Lemma 6.1. Suppose that g £ A / f°°(B n+1 , (0, T)) is a solution to the normalized Ricci- 
DeTurck flow which is e-close to h, where h is an n- admissible Einstein metric and e is 
sufficiently small. Then 

d t \g - h\ 2 < g^ViVjlg - h\ 2 - (2 - e)\V(g -h)\ 2 + (4 + e + b{n)n)\g - h\ 2 , (6.1) 

where b(n) is a constant depending on n, and where any constant c(n)e is replaced with e, 
and norms and covariant derivatives are with respect to h. 

Proof. Choose h- normal coordinates at any point p and write hij = 5ij, then diagonalize g at 
p to write gij = Xigij, where Aj > 0. Note that we drop Einstein summation convention as we 
use this expression in the remainder of the proof. 

Note that for a metric g e-close to h, the eigenvalues of g with respect to h satisfy (l+e) _1 < 
Ai < 1 + e. 

Set Z = g — h. We compute using equation (2.2) 
d t \Z\ 2 = 2^2(gu - hu)dtgii 

i 

= 2Y J {9u ~ ha) (g ab V a V bgil + 2g ab g ip h pq E iaqb 

i 

+2{ha - go) + 2g iia ab (h ab - g ab ) + [g^ 1 * g' 1 *Vg* V^J , 

where E represents the curvature deviation from hyperbolic space: E = R + R cc . 

We now estimate exactly as in [29]. The only new ingredient are the curvature error terms. 
Given an index i, we find the expression for 2^ i {ga — ha)dtgn contains 

2^22{ gii - hii)g ab g ip h pq E ia g b = 4^( 5ii - hu)g ab g iq E iaqb 

i i,q 

i 

= -4^2(X l -l)X i g ab E ailb 

i 

= -4^(A,-l)A i ^AT 1 ^ i 

i 3 
= -4^(A J -l)A^(A7 1 -l) J B jW 

i j 

where in the last line we used the fact that since h is Einstein ^ ■ Eijji = 0. Noting from 
symmetries of the curvature tensor that Ena = 0, we find 

-4^(Ai - l)Ai ^(AJ 1 - l)E m = 4^A i Aj 1 (A i - 1)(A, - l)E m 
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Estimating using the /i-norm yields 



-^(Ai-ljA^AT 1 -!)^. 



4j^A i AJ 1 (A i -l)(A i -l)E W i 



< ^max|£ , jj J j|^^ &(n)|Z| 

< r)b(n)\Z\ 2 . 



□ 



In order to proceed we must replace the McKean inequality used in [29] for the infinimum 
of the I? spectrum of the hyperbolic Laplacian on functions with its counterpart for a con- 
formally compact Einstein metric. By a result of Lee [18], if h is at least C 3 ' a conformally 
compact with smooth conformal infinity, and the conformal infinity has positive Yamabe 
invariant, then Xo(h) = n 2 /A. 

We now obtain I? control of \g(t) — h\ by modifying [29, Theorem 3.1]. Let rj(n) = ^ , 
where b(n) is the constant appearing in the previous lemma. We assume that for < r/ < r](n) 
that we have a fixed 77-admissible Einstein metric. 

Theorem 6.2. Let n > 3. There exists 5 = 5 (n) > such that if g £ M°°(Br, [0, T)) 
is a solution to the normalized Ricci-DeTurck flow with g = h on dBn(0) x [0, T) and 
su Ps fl (o)x[o,T) \g-h\<6 , then 

[ \g(t) - h\ 2 h dvol h < e~ at [ \g(0) - h\\ dvol h . 



a > 1/4. 



Proof. Assume that 5q is so small that g is e-close to h. 
We compute using the Lemma 6.1, 



d t \Z\ 2 dvol h </ s tJ ViVj|Zr-(2-e)|VZr + (4 + e + 77)|Zrdvolh 

JB R (0) Jb r (o) 

After integrating by parts and estimating as in [29], we obtain 



d t \Z\ 2 dvol h </ -(2-e)|VZ| 2 + (4 + e + 6(n)r/)|Z| 2 dvolh, 

We now use the eigenvalue estimates described above. Kato's inequality implies that 
|V|Z|| 2 < |VZ| 2 . Further, 

/ |V|Z|| 2 dvolh > X (B R (0)) [ |Z| 2 dvol h 
Jb r (o) Jb r (o) 
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where X o (B R (0)) > \ (h) = So now, 



I -(2- e)\VZ\ 2 < [ -(2-e)|V|Z|| 2 dvol h 

JBn(O) JBnfO) 



<B R (0) JB a (0) 

f n 2 
< / -(2-e) — |Z| 2 dvol h 



B R (0) 

~ + e) I |Z| 2 dvol h . 



2 



Br(0) 



Finally, we obtain 



-y +e + 4 + b(n)T] 



where F = J Br{q) |Z| 2 dvol h . 

2 

We must have —\ + e + A + b(n)r] < in order to get an exponential decay estimate. 
Thus we can take e so small that (say) e < 1/8. Since 6(71)7/ < 1/8 and n > 3, we see 
-^ + e + 4 + 6(n)r ? <-i + i = -i. ' □ 

We now have all of the new ingredients needed to prove Theorem B. Let < ry < rj(n) and 
h be an 77-admissible Einstein metric. We proceed exactly as in [29, Theorem 3.4]. 



Appendix A. More on Linear Parabolic PDE theory 

In this appendix we give more detail surrounding linear parabolic theory on conformally 
compact asymptotically hyperbolic manifolds. Our approach to understanding these operators 
is based on the edge heat calculus developed in [1]. Note that in this appendix we deal 
exclusively with the 0-case but the arguments generalize in a straightforward manner to the 
full complete edge case. 

The point of view we adopt is that for a second order uniformly degenerate elliptic operator 
with time- independent coefficients, we can explicitly construct the heat kernel as a poly ho- 
mogeneous distribution on an appropriate manifold with corners that covers M x M x R + . 
In this section we will first describe this blow up space. We then proceed to discuss the heat 
kernel as constructed in [1]. We then prove several mapping properties of these kernels. We 
conclude by proving Schauder type estimates. 

We now introduce the appropriate blow up spaces for the construction of the heat kernel. 
First we define the 0-double space: M 2 , originally introduced in [22] for the elliptic edge 
calculus. This is a manifold with corners that covers M 2 , and is obtained by introducing 
polar coordinates around the submanifold 

dM x b dM = {(w, w') £ dM x dM : w = w'}. 

So M 2 = [M x M; dM x B dM]. This will introduce three new boundary hyper surf aces; 
following Albin we denote these by B\\ (the front face), Bq\ (the right boundary) and B\$ 
(the left boundary). We denote the blowdown map 

P e : Ml -> M 2 , 

and the edge diagonal by 

diag e = ^(diag^M x B dM). 
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We describe the edge double space in terms of coordinate charts. In the interior of M 2 we 
may use the usual coordinates 

{(x,y),(x',y')) = ((,('), 

where y will always denote coordinates along B and z will always denote coordinates along 
F. We will favour the following projective coordinates for Mg, defined away from B\q and 
that express the edge diagonal easily are given by 

, \ ( x' y' -y 
{x,y,z), s := —,v := 

\ X X 

Note that in these coordinates, s = is a defining function for Bqi and x = for the front 
face (away from Bio). By reversing the roles of x and x' in the obvious manner, one may 
obtain a second chart covering the remainder of M e 2 . 

We now introduce the heat space iJM 2 . This is given by a parabolic blow up of the 
manifold Mg x ]R + along the submanifold diag e x {0}. This gives us a number of new boundary 
hypersurfaces. We keep Albin's notation for these, illustrated in Figure 1. 




We now discuss the coordinate systems we can use on _BM e 2 . In what follows we work away 
from .Bio,o (i- e - away from x = 0). Near -611,0, an d away from i?oo,i we can use 

(<-4 s ' : =^' : ~)- : =' I/2 )- 

Near -Bn,o and the 'top' of -Boo, 2 we may use 

((S,U),C',t) := {{^^)^^ 1 ' 2 ) ■ (A.2) 

In these coordinates, x' = defines the -Bn^ and r = defines -Boo, 2- Finally, near -Bn,o and 
the 'bottom' of -Boo, 2, close to -Boo,i we appear to need to introduce another coordinate system. 
However, we observe that this region is reached using the above coordinates as \ (S, U)\ — > +00. 
We will soon see that our heat kernels vanish to infinite order along this boundary. 
We will denote the full blow down map /3 : HM^ M 2 x R+. 
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Given a manifold with corners M, C°°(M) denotes functions on M that are smooth in 
the interior and smooth up to all boundary hyper surf aces. The space C°°(M) will denote 
smooth functions vanishing to all orders at the boundary hyper surf aces. If T denotes a list 
of boundary hypersurfaces then Cj?(M) denotes smooth functions vanishing to all orders at 
all boundary hypersurfaces except those in J 7 ; at the other hypersurfaces we demand the 
functions are smooth up to the boundary. 

We will also need to define sets of functions that have asymptotic expansions at the bound- 
ary hypersurfaces. Let M be a manifold with corners with boundary defining functions Xj. A 
distribution u is polyhomogeneous conormal 1 if: 



U r 



Pj 



^2a jjP (x,y)x s i (\ ogx y 



YOO p- 



where ctj tP G C°°(M). We'll denote the set of such distributions <A* h . We can also restrict the 
set of exponents that may occur above. Define an index set to be a discrete subset E C C x No 
such that 

(1) if (sj,pj) G E and \(sj,pj)\ — > oo, then Re(sj) — > oo. 

(2) if (s,p) G E then (s + k,p - I) G E for any k,l G N,Z < p. 

Given a set of index sets £ for each boundary hypersurface, we denote by A^ h the set of 
polyhomogeneous conormal functions with exponents ranging in £. Note that we will use a 
few special notations for index sets. The empty set will denote the index set for a function 
vanishing to all orders along a hypersurface. A single number n G No will denote the index 
se t {(j)0) : j G N,j > n} of functions vanishing to order n. Note that the index set {0} 
represents functions smooth up to the hypersurface. For more details about operations on 
these sets, see the concise review in [22, Appendix A]. 

A.l. The heat kernel of a uniformly degenerate elliptic operator. Let L be a second 
order uniformly degenerate elliptic operator. We consider a heat type equation 

(d t -L)u(C,t) =0 

«(C,o) =/(C), 

where / G T(M; £ ) is a smooth section of a vector bundle £. 

The heat kernel of L is a distribution on M 2 x R + so that the solution to the above problem 
is given by: 

u(C,t)= / MC,C^/(C')dvoi g (C'). 

Here h formally satisfies: 

(dt-L c )h(C,C,t) =0 

fc(C,C',o) =s((-C), 1 j 

We will see that h = (3*H, where H is a polyhomogeneous distribution on HM 2 . 

The actual construction of this distribution is done for half-densities, so that it makes sense 
to compose operators. We briefly review Albin's construction of the heat calculus. We define 



^See [11] for a discussion and to make the meaning of ~ precise. 
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a weighted bundle of half-densities D := p 00 2 2 p 11 2 Q}' 2 {HM 2 ). Kernels of operators in the 
heat calculus are elements of 

K k ' l (M, D) := p^ 2 p\ lfi C% 02>Bll {HMl-D). 
The action of a kernel Ka in K k)l on smooth half-densities is given by 

A(f)(£,t)= [ /3*K A ((,(',t)f(('). 

We'll denote the operator A acting in this manner by A £ ^e'jjeaf 
Albin proves: 

2 

Theorem A.l. If L is the scalar Laplacian of a exact edge metric, then A € ^ e Heav w ^ iere 
A is the heat operator of dt — L. 

We note that Albin's construction is closely modeled on the work of Melrose [26], and 
generalizes in a straightforward manner to general second order uniformly degenerate elliptic 
operators. Furthermore, Melrose also considers the case of elliptic operators between bundles 
[26, Theorem 7.29] with diagonal principal symbol. Thus we have 

Theorem A. 2. If L is a uniformly degenerate elliptic operator with diagonal principal symbol, 

2 

then A £ ^ ^ Heav w here A is the heat operator of dt — L. 

We now give a brief indication of the proof of theorem and refer the reader to [1, 26] 
for further detail. We work in HM 2 with the ansatz that the solution already vanishes to 
infinite order at -Bio,0) ^oi,o> an d £?oo,i- To deal with the rest of the equation and boundary 
hypersurfaces involves three main steps. First, an initial parametrix is constructed by pulling 
the heat equation back to HM 2 in coordinates near the blown up diagonal. As -Boo, 2 fibres 
over the diagonal, we find that the equation restricts to a Euclidean type heat equation on each 
fibre with smooth coefficients in the variables along the fibre. Thus we may progressively solve 
away the Taylor series at -Boo, 2 with control of the asymptotics down to £>ii,o- This handles 
the initial condition. The second step is to progressively solve away the Taylor series at B\\^ 
using the heat kernel of hyperbolic space (recall 0- metrics are asymptotically hyperbolic). 
The result of these two steps is a parametrix solving the heat equation to infinite order at all 
boundary hypersurfaces. To improve the parametrix to an actual inverse requires an argument 
involving Volterra operators and is given in [26, Proposition 7.17]. 

Finally we note that the construction above also works when the background metric is 
p oly homogeneous . 

A. 2. Mapping properties. In this section we study the action of the heat kernels in ^g'^ eaf 
above on functions, using Melrose's pushforward theorem. Figure 2 introduces some important 
notation. 

We identify functions and half-densities on M 2 x M + and the factors M x M + and M by 2 
f(x,y,x',y',t) O f{x,y,x\y\t)x- ( ^{x')~ k! ^\dxdydx'dy'dt\ l l 2 , 

(n+1) i Icy 

f(x,y,t) <-> f(x,y,t)x —\dxdydtyi\ 

(n+l) , 



n+1 



{n-t L) I/O 

f(x,y) f{x)x —\dxdy\ L/z 



2 Here we omit the smooth factor ^/det g in the densities that follow. 
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M*R + M 

Figure 2. Definition of various maps 

From [1, page 13] an element of A G ^^ eai lias an integral kernel that may be written as 

/°oo 2 2/°ii o 2 k ■ u, where k is a function that vanishes to infinite order at -Bio,o> -Boi,o> and -Boo,i, 
and is smooth up to the boundary at i?oo,2 and -Bn,o, and v is a smooth section of Q}I 2 {HM 2 ). 

2 

An operator A G ^ e 'jj eat acts on half-densities by 

(Af)(x,y,t)x-^\dxdydt\ 1 / 2 = (fo), [p^p^ku ■ (J3 R )*(ftf, \dx'dy'\ 1 ^ . 

(A.4) 

To relate these half-densities, let us work in coordinates near Bn^ and -Boo, 2- We may take 
v = \dSdUdx'dy'dT\ l l 2 . Pulling our standard half-density on M 2 x R + back we find 

p*{x-^{x')-^\ dxdydx'dy'dt] 1 / 2 ) 

= (1 + 5r)-^(rE0-^(rE0-^(2(x0 n+ V n+2 ) 1 / 2 |d5dC/dx'(iy / dr| 1 / 2 

The factor \/2(l + <SV)~( n+1 )/ 2 is smooth and uniformly bounded, so we omit it hereafter. 

In order to apply Melrose's push-forward theorem, we must work with smooth b-densities. 
Here is how to arrange this. We multiply both sides by the half density x ~ [dxdydtl 1 / 2 , 
and noting 

PKx-^ldxdydt^^^iix'j-'^ldx'dy'l 1 / 2 ) = p* (x^ {x'Y^^xdydx' dy' dt\ 1 / 2 ), 
we find the action on smooth densities is given by 

(Af)(x,y,t)x-^\dxdydt\ = (poo^PnV^ ' &d ^) ■ 

Finally we introduce a total defining function on both sides of this equation to ob- 
tain 6-densities, which we denote with a subscript b. In this case we have v 2 = 
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(pio,oPoi,oPll ) oPoo ) lPoo ) 2) 1 u 2 , and a b := ±\dxdydt\. Now 
{Af)(x,y,t)xtx-^ a b 

= {Pl)* (pio,oPn n o/ 9 oi,i/ , oo,iPoo,2^ ■ WC/O^l/)) v b ) 
We now apply this to the following 

Proposition A.3. Let A G ^% eaV If f G A* hg {M) then Af G ^ 0) (M x M+). 

Proof. First, let us consider the b-map As no boundary hypersurface is mapped to a 
corner of M, /3r is a b-fibration. It is easy to check that the exponent matrix for this map is 

Bio,o Bn,o ^01,0 A)0,i -^00,2 
ax 1 1 

As a consequence, if / G >^ fl (M) with index set F, /3* R f G ^JJ^' 0,0 ' 0} (FM e 2 ), by the 
pull-back theorem [22, Proposition A. 13]. 

The function k is poly homogeneous with respect to the index set {0, 0, 0, 0, 0}. Accounting 
for the powers of the defining functions we obtain the index set {0, —re, 0, 0, 2}. The index set 
for the product of this expression with the pull-back of / is then Q = {0, —re + F, 0, 0, 2.}. 

The map (3l is a b-fibration. The exponent matrix for this map is 









Boi,o 


-Boo.i 


B { 


Hi 





1 


1 








H 2 











1 


1 



In the above table we have the labeled hypersurfaces of M x M + in the following manner: ~Hi 
represents t = and H2 represents x = 0. We now apply Melrose's pushfoward theorem [22, 
Proposition A. 18]. Note that the integrability condition is met at -Bio.o as Re{G{B\Q$)) > 0. 
Now the index set for Hi = g(B 11)0 )Ug(B 01fi ) = -n + F and U 2 = £(£oo,i)U£(A)0,2) = 2. 
Note that H2 = 2 is not surprising since r 2 = t, from the parabolic blow up. 

The calculation here provides the index sets needed to for computing the asymptotics of 

(Af)(x,y,t)xtx-^ a b . 

Canceling the powers of the defining functions, and returning to the identification of densities 
with functions now shows that Af G A { p ^f(M x R+). □ 

We use the above proposition primarily in the form 

Corollary A.4. If f £ x^C°°{M) then Af G x^C°°(M T ). 

Proof. The only point that we have to be careful about is that the previous theorem only 
guarantees an expansion in powers of r = \ft. However we can obtain full smoothness in 
t by using the fact that Af solves the heat equation and is already smooth in the spatial 
derivatives. □ 

Corollary A. 5. If f G x^C°°{Mt) and H denotes the time convolution of the heat operator 
of A, then Hf G x^C°°(M T ). 

We conclude this section with a proposition that we will need in the finer regularity analysis 
of the Ricci flow. This shows that the commutator of a 6-vector field with an element of the 
heat calculus remains in the calculus, and thus has the same mapping properties. 

Proposition A. 6. If A £ ^ 2 ^ Heat { x ) and V b is any b-vector field, then [A,V b ] G ^° Heat (X). 
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Proof. The proof is similar to [22, Proposition 3.30], adapted to the heat calculus setting. We 
sketch the proof here. Return to the action on half-densities, equation (A. 4), and suppose 
that / is a smooth half-density vanishing to all orders at the boundary hyper surf aces. Now 
suppose for simplicity that Vf, = d y is a 6-vector field. After an integration by parts, we may 
write 

(d y Af)(x,y,t) - (Ad y f)(x,y,t) 



where dy, is the adjoint of d y > under the measure. The key now is that while each of d y and d y * 
lifts to a vector field singular near -Boo.2> their sum cancels this behaviour. Indeed, computing 
in the coordinates defined in equation (A. 2), we find that 



Pii(dyt) = —du + d y i + smooth function. 



X'T 

Consequently, fi* L (d y ) + /3|j(c^) does not affect the asymptotics of the kernel, and [A, Vj,] S 

We conclude with a discussion of the main existence theorem for the inhomogeneous Cauchy 
problem: 

~ (d t -L)u(C,t) =f(C,t) 

u(C,0) =0, (A - 5j 

where / € x^C e ' 2 (M T ) and L is a second order uniformly degenerate elliptic operator with 
coefficients in C" . 

Theorem A.7. Suppose L is a second order uniformly degenerate elliptic operator with 

a - 

time-independent coefficients. For every f E x^Ce ' 2 (My) there is a solution u to (A. 5) 
in x^C e ' 2 (M) . Moreover, u satisfies the Schauder-type estimate 

INI 2+a ,2+« <^||/|| a , s (A.6) 

x^C e ^{M T ) ^C e 2 {M T ) 

Proof. By Duhamel's principle, a solution to (A. 5) is given by 

n(C, t)= [ [ h(C, C, t - t')f((', dvol g (C')< (A.7) 

J JM 

where h is the heat kernel of the heat operator e tL , where e tL G ^^'% ea f 

We now discuss the estimates. The case for nonzero weight fi reduces to the unweighted 

2+a 2+a 

estimate, as to solve the inhomogeneous problem with u £ x^C e ' 2 (M) amounts to solving 

(d t -x-»Lx»)u'((,t) = f'((,t) 

for with v! and /' in appropriate unweighted spaces. The kernel of the conjugated operator 
x~^Lx^ has precisely the same asymptotics as the kernel of L, as may be seen by working in 
the coordinate systems (A.l) and (A. 2). So it suffices to check the mapping properties when 
H = 0. 
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The strategy is now to cut up the space HM 2 . Consider a function <fi equal to one in a 
tubular neighbourhood of i?oo,2 and vanishing outside a slightly larger tubular neighbourhood. 
We may write the heat kernel as 

h = hi + h 2 :=<f>h + (l- <j>)h. 

To prove (3.2), it will suffice to estimate both 

«i(C,*)= / / ^(CC'^-O/CCOdvoi^cO^' 

J o Jri 

and ^ 

u 2 (C,t)= [ [ h 2 ((X',t-t')f((',t')dvo\ g ((')dt' 
Jo Jm 

Regarding the estimate for u 2 , we view h 2 as a poly homogeneous distribution vanishing 
to infinite order at -Bio,0) -E?oi,o 5 -£?oo,i) -E?oo,2 anc i smooth up to -Bii,o- The estimates are then 
checked in each of the coordinate systems (A.l) and (A. 2) in a lengthy but straightforward 
manner. For example, under the coordinate change (A. 2), the 0- vector field xd x lifts to 
(t" 1 + S)ds, consequently 

(xdx) 2 u(x,y,t) = II (xdx) 2 h 2 ((,t',t-t')f(t',t')dvo\ g ((')dt' 
Jo Jm 

= Jo I (( T " 1 + 5 )^) 2 ^(5,^x / ,y / ,t-r)/( 2 ; / ,2/ / ,r)dvol g (x , ,y / )^- 

We can now estimate the L°° norm and Holder seminorm of (xdx) 2 u(x, y, t) using the fact 
that h 2 is smooth and vanishes to infinite order in r to absorb the apparent singular factor 
of r. The same is true for the other derivatives. We omit the precise estimation but refer the 
reader to [3] for similar estimation. 

Regarding the estimate for u%, recall the Whitney decomposition outlined in Section 3.1. 
Let tpi : Bt — > (-Bj)t be an affine map taking a 'standard' parabolic cylinder over the ball 
with centre (1,0) and radius 1 in the right half plane model of hyperbolic space to the ball 
with centre (xi,yi) and radius ^xf. 

tpi(v,w,t) = (xiV,yi + XiW,t). 

Consider the pullback of u\ under tpf. 

(ip*u)(v,w,t) =u 1 (x i v,y i + x i w,t) 

-111- 

-a- 

These charts essentially provide the bookkeeping for rescaling of the equation to a fixed 
parabolic cylinder in M n+1 x [0, oo). Further, uniformly degenerate vector fields pull back under 



/ XjV - x' yj + XjW - y' , , _ A 

\ X '{t-t')y^ x'it-f) 1 / 2 ,x,y ' ) 

■f{x',y',t'){x>)- n - l &vo\z{x',y')dt' 
( v — v w — w „ „ A 

[ v(t-tr^ v(t-t>yi^ m + XiW > 1 - 1 ) 

■ f(xiV, yi + X{W, t')v~ n ~ l dvol g (?;, w)dt' 
( v — v w — w _ _ A 

• {ip*f){v,w,t')v- n ~ 1 dvol g {v, w)dt' 
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ipi to differential operators with uniformly bounded coefficients on the standard cylinder, and 
in particular a uniformly degenerate parabolic operator is pulled back to a uniformly parabolic 
operator. This reduces the estimate to the 'classical' parabolic case, which may be found in 
[17]. 

□ 
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